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Abstract —The use of optical spatial modulation (OSM), which 
has been recently emerged as a power and bandwidth efficient 
pulsed modulation technique for indoor optical wireless com¬ 
munication, is proposed as a simple, low-complexity means of 
achieving spatial diversity in coherent free space optical (FSO) 
communication systems. In doing so, this paper makes several 
novel contributions as follows. It presents a generic analytical 
framework for obtaining the Average Bit Error Probability 
(ABEP) of uncoded OSM with coherent detection in the presence 
of turbulence-induced fading. Although the framework is general 
enough to accommodate any type of models based on turbulence 
scattering, the focus in this paper is the H-K distribution. 
Although this distribution represents a very general scattering 
model valid over a wide range of atmospheric conditions, it is has 
not been considered in the past in conjunction with FSO systems 
possibly because of its mathematical complexity. The proposed 
analytical framework yields exact performance evaluation results 
for MIMO systems with two transmit and an arbitrary number of 
receive apertures. In addition, tight upper bounds are derived for 
the error probability for OSM systems with an arbitrary number 
of transmit apertures as well as for convolutionally encoded 
signals. The performance of OSM is compared to that of well 
established coherent FSO schemes, employing spatial diversity 
at the transmitter or the receiver only. Specifically, it is shown 
that OSM can offer comparable performance with conventional 
coherent FSO schemes while outperforming the latter in terms of 
spectral efficiency and hardware complexity. Various numerical 
performance evaluation results are also presented and compared 
with equivalent results obtained by Monte Carlo simulations 
which verify the accuracy of the derived analytical expressions. 

Index Terms —average bit error probability, atmospheric tur¬ 
bulence, coherent detection, free space optical communica¬ 
tion systems, H-K distribution, multiple-input multiple-output 
(MIMO) systems, optical spatial modulation. 

I. Introduction 

Free-space optical (FSO) communication systems have re¬ 
cently attracted great attention within the research community 
as well as for commercial use. FSO systems can provide ultra- 
high data rates (at the order of multiple gigabits per second), 
immunity to electromagnetic interference, excellent security 
and large unlicensed bandwidth i.e. hundred and thousand 
times higher than radio-frequency (RF) systems, along with 
low installation and operational cost 0. 
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The challenge in employing such systems is that FSO links 
are highly vulnerable due to the detrimental effects of atten¬ 
uation under adverse weather conditions (e.g. fog), pointing 
errors and atmospheric turbulence. One method to improve 
the reliability of the FSO link is to employ spatial diversity, 
i.e. multiple-lasers and multiple-apertures to create a multiple- 
input multiple-output (MIMO) optical channel. Because of its 
low complexity, spatial diversity is a particularly attractive 
fading mitigation technique and performance enhancements 
have been extensively studied in many past research works in 
the field of FSO communications 0-0. 

In order to evaluate the impact of atmospheric turbulence on 
the performance of OSM, accurate models for the fading dis¬ 
tribution are necessary. For example the lognormal distribution 
is often used to model weak turbulence conditions whereas the 
negative exponential and the K-distribution are used to model 
strong turbulence conditions 0. Other more general statistical 
models have also been proposed to model scintillation over 
all turbulence conditions, including the Gamma-Gamma 0, 
the lognormal-Rice (or Beckmann) 0 the homodyned K 
distribution (H-K) ffTOl and the I-K llTTl - liT3l distributions. All 
these three models are based on the argument that scintillation 
is a doubly stochastic random process modelling both small 
and large scale turbulence effects. Besides, they agree well 
with measurement data and simulations for a wide range of 
turbulence conditions. 

In this paper, the H-K distribution is adopted to model 
turbulence-induced fading. The main reason for this choice 
is the fact that this distribution is based on a very general 
scattering model which is valid for a wide range of atmo¬ 
spheric conditions. It is also noted that the H-K distribution 
generalizes existing models such as the K-distribution. The H- 
K distribution models the field of the optical wave as the sum 
of a deterministic component and a random component, the 
intensity of which follows the Rice (Nakagami-n) distribution. 
The average intensity of the random portion of the field 
is treated as a fluctuating quantity im. It is important to 
underline that, to the best of our knowledge, in the open tech¬ 
nical literature there have been no papers published analyzing 
and evaluating the performance of FSO systems over such 
channels, because of the complicated mathematical form of 
their respective probability density functions (PDF). 

Depending on their detection, FSO systems can be classified 
into two main categories, namely coherent (heterodyne detec¬ 
tion) and non-coherent (direct detection) systems. Coherent 
FSO systems have the information bits encoded directly onto 
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TABLE I 

List of Mathematical Notations 

j 2 = — 1 denotes the imaginary unit 
\z\ denotes the magnitude of the complex number z 
denotes the real part of the complex number z 
denotes the real part of the complex number z 
f{x) = o[g{x)] as x xq if = 0 

|| • ||^ denotes the square Frobenius norm 
(•) T denotes the matrix transpose 

* denotes convolution 
E(-) denotes expectation 

fx (•) denotes the Probability Density Function (PDF) of the random variable X 
Fx(-) denotes the Cumulative Distribution Function (CDF) of the random variable X 
Aix(‘) denotes the Moment Generating Function (MGF) of the random variable X 
I a (•) is the modified Bessel function of the first kind and order a CD eq. (8.431)] 

K a (• ) is the modified Bessel function of the second kind and order a (T| eq. (8.432)] 
r (x) = exp(— t)t x ~ J df is the Gamma function (T] eq. (8.310/1)] 

Q{x) = —exp(— t 2 /2)dt is the Gauss Q-function 
Wp.q(-) is the Whittaker function [T] eq. (9.220)] 

Pr{ } denotes the probability operator 

• denotes estimated value at the receiver side 


the electric field of the optical beam. At the receiver, a local 
oscillator (LO) is employed to extract the information encoded 
on the optical carrier electric field. On the one hand, coherent 
FSO systems can provide significant performance enhance¬ 
ments due to spatial temporal selectivity and heterodyne gain 
in comparison to direct detection systems. Moreover, they are 
more versatile as any kind of amplitude, frequency, or phase 
modulation can be employed. On the other hand, coherent 
receivers are more difficult to implement as the LO field should 
be spatially and temporally coherent with the received field. 

Recently, the so-called optical spatial modulation (OSM) 
has emerged as a power- and bandwidth-efficient single-carrier 
transmission technique for optical wireless communication 
systems m-sm. This spatial diversity scheme, initially 
proposed in fl7l and further investigated in IT8l . 1 1 9| , employs 
a simple modulation mechanism that foresees to activate just 
one out of several MIMO transmitters at any time instant and 
to use the index of the activated transmitter as an additional di¬ 
mension for conveying implicit information. It has been shown 
that OSM can increase the data rate by base two logarithm of 
the number of transmit units M- Also, OSM can increase the 
data rate by by a factor of 2 and 4, respectively, as compared 
to on-off keying (OOK) and pulse position modulation (PPM) 
urn CD]. It is underlined that such performance gains are 
obtained with a significant reduction in receiver complexity 
and system design. 

Because of the above mentioned advantages of OSM over 
other more conventional transmission schemes and given the 
wide applicability of FSO, it is of interest to investigate the 
potential performance enhancements obtained by incorporating 
OSM in FSO systems. However, in general this research 


topic has not been dealt within our research community. 
Only recently, there have been papers published in the open 
technical literature dealing with performance analysis studies 
of FSO systems employing spatial modulation and operating 
in the presence of atmospheric turbulence, e.g. see fl20l 
and ED. Specifically, in 120), the combination of subcarrier 
intensity modulation and spatial modulation with receiver 
diversity was proposed to enhance the performance of intensity 
modulated direct detection (IM/DD) FSO systems. In ED . 
another IM/DD based system FSO system which combines 
antenna shift keying with joint pulse position and amplitude 
modulations was considered. For this system, which was 
denoted as spatial pulse position and amplitude modulation 
(SPPAM), the atmospheric turbulence channel was modeled 
as log-normal or Gamma-Gamma distributions and was eval¬ 
uated, in terms of bounds, for uncoded and coded signals. 
ABEP performance evaluation results have shown that SPPAM 
offers a compromise between spectral and power efficiencies 
as well as a certain degree of robustness against atmospheric 
turbulence. Despite these two papers which deal with non¬ 
coherent detection schemes, the potential enhancements of 
OSM on the performance of FSO systems with coherent 
detection still remains an open research topic which, to the 
best of our knowledge, has not been addressed so far in the 
open technical literature. 

Motivated by the above, in this paper we present for the first 
time a generic analytical framework which can be used to ac¬ 
curately obtain the performance of outdoor OSM with coherent 
detection in the presence of turbulence-induced fading. More 
specifically and within this novel analytical framework, the 
main novel research contributions of the paper are as follows: 
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• New analytical expressions for the ABEP of coherent 
OSM under turbulence conditions modeled by the H-K 
distribution are derived. When the transmitter is equipped 
with two apertures the resulting analytical expressions 
are exact, whereas for an arbitrary number of transmit 
apertures tight upperbounds can be obtained. 

• Error probability performance bounds for coded OSM 
systems are derived and the performance enhancements 
when channel coding is employed are presented and 
analyzed. 

The error probability performance of OSM is also compared 
to that of conventional FSO schemes with transmit or re¬ 
ceive diversity only, i.e. when Maximal Ratio Combining 
(MRC), Selection Combining (SC) or Alamouti-type Space- 
Time Block Codes (STBC) are employed. It is noted that 
the theoretical analysis is substantiated by comparing the 
theoretical and equivalent simulated performance evaluation 
results obtained by means of Monte Carlo techniques. 

The paper is organized as follows. After this introduc¬ 
tion, Section [Tl] outlines the system and channel models. In 
Section [III] analytical expressions for the ABEP of uncoded 
OSM systems are presented. Asymptotic ABEP expressions 
are also derived, wherefrom the diversity gain of coherent 
OSM can be readily deduced. The performance of coded OSM 
systems is discussed in Section [Tv] In Section [V] the various 
performance evaluation results and their interpretations as well 
as comparisons are presented. Finally, concluding remarks can 
be found in Section [VI] Notations: A comprehensive list of 
all mathematical notations used in this paper can be found in 
Table U 

II. System and Channel Model 

In this section, a detailed description of the OSM FSO sys¬ 
tem model , i.e. transmitter, channel and receiver is provided. 
Moreover, the H-K distribution is introduced and analytical 
expressions for its parameters in terms of equivalent phys¬ 
ical parameters of the turbulence phenomenon, such as the 
refractive-index structure parameter, optical wave number, and 
propagation path length, are derived. 

A. Preliminaries 

Let us consider a M x N MIMO FSO system with M 
transmit units (lasers) and N coherent receivers. It is assumed 
that the receiving apertures are separated by more than a 
coherence wavelength to ensure the independency of fading 
channels. The basic principle of OSM modulation is as follows 

m, G2): 

i) The transmitter encodes blocks of log 2 (M) data bits into 
the index of a single transmit unit. Such a block of bits is 
hereafter referred to as “message” and is denoted by b m , 
Vm = 1,2It is assumed that the M messages are 
transmitted with equal probability by the encoder and that the 
related transmitted signal is denoted by E m = E m exp(j^f, m ). 
During each time slot, only one transmitter £, where l = 
1,2, ...,M is active for data transmission. The information 
bits are modulated on the electric field of an optical signal 
beam through an external modulator. During this particular 


time slot, the remaining transmit lasers are kept silent, i.e. 
they do not transmit. 

ii) At the receiver, the incoming optical field is mixed 
with a local oscillator (LO) field and the combined wave is 
first converted by the photodetector to an electrical one. A 
bandpass filter is then employed to extract the intermediate 
frequency (IF) component of the total output current. Finally, 
a A r -hypothesis detection problem is solved to retrieve the 
active transmit unit index, which results in the estimation of 
the unique sequence of bits emitted by the transmitter. 

B. Receiver Structure 

The received electric field at the aperture plane of the ?i-th 
receiver after mixing with a LO beam, can be expressed as 

ED, m 

^n(^) = y/ < ^Pt^()E rn hm,n COs{idQt ~\~ (f>rn,n H - 4*bm ) ^ ^ 

+ V ZPloZo COS (uJLOt)- 

In the above equation, P* is the transmit laser power, Zq is the 
free space impedance, h m ,n and <j>m,n denote the magnitude 
and the phase of the complex channel coefficient between 
the m-th transmit and the 77 -th receive aperture, respectively. 
Furthermore, Plo denotes the power of the local oscillator, 
WiO = wo + ijJif where tu 0 and uijp are the carrier and the 
intermediate radian frequencies, respectively. 

The output current of the n-th photodetector can be math¬ 
ematically expressed as If 22 ) . |23| 

in(t) = ■!■[«5„(i)] 2 (2) 

where R = t]q e /(huo) is the responsivity of the photodetector 
with q e = 1.6 x 10 -19 Cb is the charge of an electron, 
h = 6.6 x 10~ 34 J • s is the Planck constant, 77 is the pho¬ 
todetector efficiency, and uq = u}q/(2tt) is the optical center 
frequency. Expanding C3 and ignoring the double-frequency 
terms that are filtered out by the bandpass filter, the resulting 
photocurrent can be expressed as 

in(t) = RPtEmhm,n + RPlO 

T £ 11 \/ l 1 ; 1^1.() b'rrijl rri .ri COS^UJjpt ) (3) 

— *DC (t) + * AC CO- 

In @. * dc (0 — RPtEmh'm,n + PPlo is the DC component 
generated by the signal and local oscillator fields, respectively, 
*ac(0 = 2 Ry/P t P L o cos (ui F t-(t>m,n-4>b m ) is the AC com¬ 
ponent in the received photocurrent which, unlike for direct 
detection, contains information about the frequency and phase 
of the received signal. It is assumed that for coherent detection 
the intermediate frequency ojj p is nonzero, so that the signal 
power can be expressed as P s = 2R 2 PtPpoEm^m n 

As in l22ll - ll25ll . we also consider that Ppo Ps and thus, 
the DC photocurrent can be approximated as inc(t) ~ RPlo- 
The photodetection process is impaired by shot noise with 
variance L = 2 q e RPpoB e where B e is the electrical 
bandwidth of the photodetector. It is also noted that because 
of the large value of RPlo the photocurrent due to thermal 
noise and the dark current can be ignored li 22 l . 
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Following |f23ll and f24l . the sufficient statistics at the n-th 
coherent receiver can be expressed as 

Vn — 6Xp[j(^) m n -f- T Z n (4) 

where p = RP t /(q e B e ) is the average signal-to-noise ratio 
(SNR) and z n is the noise at the n-th receiver. Assuming that 
the LO power is large and the receiver noise is dominated by 
LO related noise terms, the Additive White Gaussian Noise 
(AWGN) model can be employed as an accurate approxima¬ 
tion of the Poisson photon-counting detection model ( 22l . (23|. 
Thus, z n can be modeled as a zero-mean unit variance complex 
Gaussian random variable ll23l . 

Similar to ll26l , it is assumed that the receiver has knowl¬ 
edge of the actual fading gains and that the total fading remains 
constant over one bit interval and changes from one interval to 
another in an independent manner. At the receiver, the optimal 
spatial modulation detector estimates the active transmitter 
index, l, at a given time slot according to m 

i = argmaxp y (y|x, H) 

( (5) 

= argmin {^/Jl || hexe \\ 2 F -2 ( y T h e xe )} 

t 

where 

- x is an M-dimensional vector with elements correspond¬ 
ing to the electrical field E m exp(j 4>b m ) that is transmit¬ 
ted over the optical MIMO channel; 

- H(t) is an N x M optical MIMO channel defined as 

H(t) = [hi,h 2 ,...,h M ] 

hn(t) exp(j^>n) ... h 1M {t) exp(j</>i M ) 

a /i2i(f)exp(.^2i) ■■■ (i2M(f)exp(j</>2M) 

h N i(t) exp(j(/)jvi) ... h NM (t) exp(j^jvM) 

(6) 

- z is the iV-dimensional noise vector; 

- p y (y|x, H) is the PDF of y conditioned on the transmit¬ 
ted vector x and the channel H; 

C. Channel Model 


random fluctuations in the turbulence parameters is modeled 
by allowing random variations in the parameter bij of the 
Rayleigh component. Following El, it is further assumed that 
follows a gamma distribution with PDF given by 


K(b) 



fray -1 
F (Olij) 


exp 



( 8 ) 


where a is the shaping parameter and represent the effective 
number of scatters and b<j t = E {bij}. Then, the PDF of the 
irradiance 1^ = \hij(t)\ 2 , . (I), can be expressed as ifTTl Eq. 

( 8 )] 


/ly(I) = 


ij/bOij) 


T (aij) 


exp - 


u 0 a 



which is actually the integral representation of the H-K dis¬ 
tribution m. It is noted that /[.. (I) cannot, in general, be 
expressed in closed form, with the exception of the special 
cases = 0 or a = 1. Specifically, for Ajj = 0 © reduces 
to the K-distribution whereas for a = 1, © reduces to a 
special case of the I-K distribution liTOl Eq. (10)]. 

The n-th normalized moment of I,j is given by GO] Eq. 
(22)] as 

E {I ij} _ A. ^ fv\ r {aij + u-k) (aijPij) 1 ' 

maV ~ + Pij) v h U7 r(ay) a 

( 10 ) 

where pij = j4?-/6oy I s the coherence parameter, defined as 
the power ratio of mean intensities of the constant-amplitude 
component and random component of the field in © IfTTi . 
E2- Using 03, the scintillation index can be readily calcu¬ 
lated as 

2 a } _ -I _ a ij ~1~ ‘^aijPij + 2 

Iij E{^'} 2 " ^{l + pijf ■ 

Under the assumption of spherical wave propagation, of., 
can be directly related to atmospheric conditions as E2 Eq. 
(7), Eq. (9)] 


A discrete scattering model is considered, where the radia¬ 
tion field of an optical wave at a particular point is assumed 
to be composed of a number of scattered components that 
have traveled different paths. Under the Ricean assumption 
CD, the complex channel path gains hij(t) between the i- 
th transmitter and the j-th photodetector can be expressed as 
hij(t) = hij(t) exp(jojt) where uj is the radian frequency of 
the optical wave and 

hij(t) = ®{hij(t)} + j $s{hij(t)} 

= A^ exp \j6ij(t)] + Rij(t) exp[j<Fy (t)] 

where the term A i3 exp (jdij(t)) is a deterministic component 
and Rij(t) exp(j$>ij(t)) is a circular complex Gaussian ran¬ 
dom variable. Hence, the amplitude Rij is Rayleigh distributed 
with parameter af ? = bij /2 IT ll Eq. (13)] and the phase <1 \j 
is uniformly distributed over [0, 2-k). Under the assumption 
of a doubly stochastic scintillation model CD, the effect of 




0.41a«2(l + 0.5tr?) ) a i < 1 
1 + 2.8/o-f 75 , o-i > 1 


( 12 ) 


where a? = 1.23Cf /c 7,/6 T.f 76 is the Rytov variance, k = 
27r/A is the optical wave number with A being the wave¬ 
length, Lij is the link distance and C nij denotes the index of 
refraction structure parameter. For FSO links near the ground. 
Cl « 1.7 x 10 -14 m -2 / 3 and 8.4 x 10- 15 m" 2 / 3 for the 

rl T-3 _ 

daytime and night, respectively B281 . Moreover, cr i « 1 and 
a i>l correspond to weak and strong turbulence conditions, 
respectively. 

Using (fl2l >. the parameters of the H-K distribution, a 
and p, can be directly related to physical parameters of the 
turbulence by following a similar line of arguments as in lfT2l . 
where similar results were derived for the I-K distribution. 
In particular, on the one hand, weak turbulence conditions 
are characterized in the H-K distribution by large values of 
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Pij. In this case the scintillation index given by (fTTb can be 
approximated as 

9 2 

erf « —, with » 1. (13) 

Pij 

On the other hand, assuming strong turbulence conditions 
where pij tends to zero, (fill can be approximated as 

9 2 

erf.. ~ 1 H-, with p^ < 1. (14) 

1J 

By comparing (fT3l> and ( fTTb with the first and second branches 
of (IT 2 b . respectively, a, :) and p 13 can be obtained as 

= 0-71 
4.88 

Plj ~ ap (1 + 0.2 al.Y 

\ l-zj f 

To the best of our knowledge, the relationship of a i3 and p^ 
with u l4 . given by (ITT and (fl6b is a novel result. 

III. Performance Analysis of Uncoded OSM 
In this section, by employing the well-known MGF-based 
approach for the performance analysis of digital commu¬ 
nications over fading channels ll29l . analytical expressions 
for the ABEP of uncoded OSM systems will be derived. 
Expressions for the diversity and coding gains of OSM systems 
are also presented, thus providing useful insight as to how 
these parameters affect the overall system performance. 


(15) 

(16) 


A. Preliminaries 


For M = 2, the conditional bit error probability (BEP) of 
OSM systems when no turbulence induced fading is consid¬ 
ered can be obtained in closed form as lf26l 

Pe( hr, h 2 ) = Q ^\\ h x - h 2 ||^ . (17) 

The squared Frobenius norm in (flTt can be expressed as 

N 

|| hi - h 2 ||2 = £ \h hn - /i 2 ,„| 2 (18) 

71=0 


where /+ n is the n-th element of 1+, Vi £ {1,2}. When 
M > 2 transmitters are considered, a tight upper bound for 
the conditional BEP of the above system can be obtained as 
M Eq. (7)] 


Pe( H ) < 

M 


M- 1 

log 2 (M) 

M 


*E E Nb(rrii, 77i2)PEP(mi —j m 2 ) 

771i = 1 7712^7711 = 1 


09) 


where PEP(toi -+ m 2 ) denotes the pairwise error probability 
(PEP) related to the pair of transmitters m\ and m 2 , where 
mi and m 2 £ 1,2,..., M, and 7V}(mi, m 2 ) is the number of 
bit which have occurred when the receiver decides incorrectly 
that m 2 instead of mi has been active. The PEP(mi —> m 2 ) 
can be evaluated as M Eq. ( 8 )] 

PEP (mi -+ m 2 ) = Q ^|| h mi - h „ l2 ||^ . (20) 


B. MGF-Based Approach 

When atmospheric turbulence is taken into account, the 
conditional error probabilities in G3 and dm need to be 
averaged over the elements of the channel matrix H in order to 
evaluate the ABEP. Without loss of generality, let us consider 
the case of a 2 x N MIMO system. Since are complex 
Gaussian random variables, the difference A„ = hi t „ — h 2 , n 
is a complex Gaussian random variable having mean equal to 
the difference of the means of and variance equal to the 
sum of variances of (+„. In order to deduce a closed form 
expression for the ABEP, it is further assumed that h^ n have 
uncorrelated real and imaginary components with the same 
variance a ^ = b n /2. It is noted that such an assumption is 
justified for link distances of the order of km and for aperture 
separation distances of the order of cm ESI, G2- For example, 
in m it was reported that for a link distance of 1.5 km, a 
wavelength of 1550 nm, an aperture diameter of 1 mm and 
photodetectors separated by as little as 35 mm, which validates 
the independence assumption. 

Consequently, A n has uncorrelated components too and its 
squared envelope, | A„| 2 , is characterized by a non-central chi- 
square PDF as follows 


f\A n \ 2 i x \bn) = ^“ eX P 



( 21 ) 


where A n = \A 2 t ne j92,n — Ai in e J01 - n \. Assuming that b n 
follows a gamma distribution with parameters a n and 
the unconditional PDF of |A „| 2 is obtained by averaging (l2ll) 
with respect to b n , i.e. 


f\A n \ 2 { x ) 


(Ctn/&0,n) an 

2r(a n ) 




To 



( 22 ) 


As was pointed out in ED, the integral in (f22l> cannot be 
solved in closed form. Nevertheless, for the special case of 
a n = 1 , i.e. when one scatterer per branch is considered, and 
by employing ED Eq. (10)], this integral can be evaluated in 
closed form as 

-^Kq (v^XEE) To (\/2-r/&o,n) , x < A 2 n 

b^-To ^2A„/6 0i T K 0 (V2x/b 0 , n ) , x > A 2 . 

(23) 

Moreover, for the special case where hi tU and h 2 , n have 
identical mean value, i.e. when A n = 0, d22l> yields the well 
known K-distribution with PDF given by 


f\A n \ 2 (x) = 


f\A n \ 2 {x) 


2 (l-a„)/2 p( an ) 




(24) 


By employing the MGF-based approach for the performance 
analysis of digital communications over fading channels, the 
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average PEP (APEP) can be obtained as 


APEP = 


1 

7T 



N r 


n 


M 


|A n | 2 



de. 


(25) 


Moreover, using the tight approximation for the Gaus¬ 
sian Q-function presented in f32l Eq. (14)] (i.e., Q(x ) ~ 
1/12 exp(—x 2 ) +1/4 exp(—2 x 2 /3)), an expression accurately 
approximating APEP can be deduced as 


1 N 

ApEp .-n[x M .(f) 

n =1 


\ IT [-^lAnl 2 



(26) 

In the following analysis, analytical expressions for the MGF 
of |A„| 2 will be deduced. Specifically, the following result 
holds: 


Proposition 1. An integral representation for the MGF of 
|A n ] 2 can be deduced as 


as 1341 


Al|A„p(s) 


(^n/^0,n) 

r(a„) 


J 

xVw,- 3 


ti 


/ J J 

j =0 


‘ItjS -f 1 


exp - 


A n s a n tj 

2 tjS + 1 6o, ra 


(31) 


where J is the number of integration points, tj are the 
abscissas and Wj the corresponding weights. In |[35l eqs. (22) 
and (23)], tj and Wj are defined as 


tj = tan 


7r 

— cos 
4 


2j ~ 1 . 
2 J ‘ 


4 J cos 2 [f cos(%i7r) + f] 


(32a) 

(32b) 


For the special case of A n = 0, it can be shown that (l27l > can 
be evaluated in closed form. Specifically, the following result 
holds: 


A4|a „| 2 (s) = 


(ttn/^0,n) 

r(a„) 


r°° uan —l 


■exp - 


A n s a n b 


(27) 


0 2 bs + 1 l 2 bs + 1 b { 


d b. 


0,n 


Proof: By employing the definition of the MGF, 
A4 |a„| 2 (s) can be obtained as 


n OO 

A^|a„| 2 (s)= / exp(-sa;)/jA„| 2 (a;)da 

Jo 


r r ^ o:+a 2 

f / exp —sx — 


(ctn/^0,n) 

2r(a n ) 7 0 


bo,r, 


2b 


X In 


A n y/x I 2 




(28) 


By changing the order of integration, the above equation can 
be expressed as 


A4|a„| 2 (s) = 


(ctn/^0,n) 

2r(a n ) 


Ot-n p oo 


r«r.-2 


exp - 


CXnb 

bo,n 


Al 


exp — sx — 


2b 


Iq 


A n ^fx 


dx 


(29) 


d b. 


The inner integral, i.e. with respect to x can be evaluated by 
employing [[33] Eq. (3.15.2.2)] as 




exp — sx — 


2b 


2b 


2sb + 1 


exp 


2A n b(2sb + 1) 



(30) 


Substituting into < [29t and after some straightforward 
manipulations, ([27} is readily deduced thus completing the 
mathematical proof. ■ 

The integral in ( 127} can be accurately approximated by 
employing a Gauss-Chebyshev Quadrature (GCQ) technique 


Corrolary 1 .For the special case of A = 0 the MGF of \ A„ | 2 
can be deduced in closed form as 


•AT|a„| 2 (s) = 


QLr\ 


2sb { 


0,n 


exp 


OL r , 


4s6i 


'0,n 


X W 




y 2sho,^ 


(33) 


This result can be readily deduced by employing the in¬ 
tegral representation of the Whittaker IF-function given in 
m Eq- (9.222)]. Moreover it is worth pointing out that ([33} 
is in agreement with a previously known result, namely the 
analytical expression for the MGF of the K-distribution. Il36l 
Eq. (4)]. 


C. Analysis of the Diversity Gain 

The diversity gain of the considered OSM MIMO system 
can be obtained by using the approach presented in j37j. 
In particular, a generic analytical expression, which becomes 
asymptotically tight at high SNR values, will be derived for 
the APEP appearing in (l25l >. as follows: 


Proposition 2. For high SNR values, (125b can be approxi¬ 
mated by 


APEP 




2tv-tr (7V+ i) 

v^r (n + 1) 



(34) 


where 

(35) 


Proof: According to |[37l Proposition 3], the asymptotic 
error performance of the OSM system depends on the be¬ 
havior of A4|A„p(s), as s — > oo. To determine an analytical 
asymptotic expression for APEP a Taylor series expansion is 
employed to approximate A4 |a„| 2 (s) as 

|A4|a„| 2 (s)| = C n \s\~ d " + o(|s| -dn ), 5 ^ oo (36) 
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where c n and d. n are parameters that determine the diver¬ 
sity and coding gains of the n-th diversity branch, respec¬ 
tively. Observe that since As/(2sb+l ) ft A/(2b) and 
l/(2sb + 1) S ft°° 1/(26s). (|27) yields 


A1|a„| 2 (s) 


(ctn/^0,n) 

2sr(a n ) 


f 


( A n cy n b \ 

ex P — —— t— I do. 


(37) 


2b bi 


0,n 


By employing (33] Eq. (2.2.2.1)], (l37l > can be solved in closed 
form yielding 


X|A„| 2 (s) 


(c^n/bo^n) 

sT(a n ) 


x K a -i 


1 2/La, 


(38) 


DO ,n 


By comparing (f38l) and (l36l > it is readily deduced that d n = 
1 and c n is given by (l35l >. Thus, by substituting ( 13 fil l into 
(l25l >. the asymptotic PEP expression can be obtained as in 
(f34l> which concludes the proof. ■ 

From (f34l i it is clear that the diversity gain achieved by 
the considered system is equal to N. It is also evident that 
the diversity gain depends only on the number of the receive 
apertures and is independent of the fading severity. This 
finding is in agreement with relevant findings reported in 
1261 and |38l , for the case of radio-frequency MIMO wireless 
systems. 

It is noted that for the special case A n = 0, i.e. when |A n | 2 
follows the K-distribution, by employing the asymptotic result 
K t (x ) ft (T(f)/2) (2/x) t l34l . c n can be further simplified 
as 

C n = -7T- (39) 

2b 0 ,n(a n - 1 ) 


IV. Performance Analysis of Coded OSM over 

TURBULENCE CHANNELS 


When coded OSM is employed, the input signal s (t) is 
first encoded by a convolutional encoder. The encoded data 
are interleaved by a random block interleaver and transmitted 
through the optical wireless channels using spatial modulation. 
It is also assumed that perfect interleaving at the transmitter 
and de-interleaving at the receiver is used. Assuming max¬ 
imum likelihood soft decision decoding, the log likelihood 
ratios (LLRs) for the i-th constellation bit when the /-th 
transmitting antenna is active are computed as Ifl4l Eq. (6)] 


LLR = log 
= log 


Pr{r = i|y| 
Pr{L = 0|y} 
Eielj e x P (-11 y 

exp (-H y 


V/ II 2 /No) 

M || 2 /Ao) 


(40) 


where £ £ {1 : M} is the set of spatial constellation points, 
£\ and £' {} are subsets from £ containing the transmitter 
indices having ”1” and ”0” at the i -th bit, respectively. The 
resulting data are finally decoded by a Viterbi decoder. 


A union bound on the ABEP of a coded communication 
system can be evaluated as ll29l 

Pub < - V P(X) V q(X, X')PEP(X, X') (41) 

n z 

X X^X' 


where P(X) is the probability that the coded sequence X is 
transmitted, g(X, X') is the number of information bit errors 
in choosing another coded sequence X' instead of X n is the 
number of information bits per transmission and PEP(X, X') 
is the pairwise error probability, i.e the probability of selecting 
X' when X was actually transmitted. 

By employing |f29l p. 510], (l-ITl) can be efficiently evaluated 
as 


i _, / ,7r / 2 

P,b<-]TP(X) / 


dN 


T[D(0),N} 


IV=1. 


(42) 


where T[D(9),N] is the transfer function of the employed 
convolutional code, N is an indicator variable taking into 
account the number of the erroneous bits and D(9) depends 
on the underlying PEP expression. Furthermore, assuming that 
uniform error probability (UEP) codes are considered and 
taking into account the symmetry property this code family 
exhibits, thus making the distance structure of a UEP code 
independent of the transmitted sequence, (l42l) can be further 
simplified as ll29ll 


Pub < 


,77/2 

ri d 

- 

L 


N= 1- 


(43) 


For M = 2, using (ITTI i. (ITSl) and Craig’s formula for the Gauss 
Q-function, i.e. Q(x) = 1/n f^ /<2 exp(— x 2 /2 sin 2 6)d0, D(0) 
can be expressed as 


where can be obtained from (l27l) . When M > 2, by 

employing IH Eq. (13)], and using a similar line of arguments 
as in the case of M = 2, D(9) can be written as 


M M 

n n M \A mi , m2 \* 

m 1 = 1 777.27^^1 = 1 



(45) 


where |A mijm2 | 2 =|| h mi - h m2 || 2 . The last MGF can be 
easily computed analytically with the help of ( 1771 ). 


V. Performance Evaluation Results and 
Discussion 

In this section the various performance evaluation results 
which have been obtained by numerically evaluating the 
mathematical expressions presented in Sections [III] and [IV] 
for uncoded and coded OSM systems operating over H- 
K turbulent channels will be presented. In particular, for 
uncoded OSM systems the following performance evaluation 
results have been obtained: i) ABEP vs. SNR for 2 x N r 
OSM systems (obtained using (l26l ) with ( 1271 ). and ( l34l > - see 
Figs. Qj [H and [3); ii) ABEP vs. SNR for 2 x N MIMO 
OSM systems, 2 x N MIMO (obtained using ( l26l ) with (ITT])). 
For the uncoded schemes, in order to validate the accuracy 
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Fig. 1. ABEP of uncoded OSM for 2 X N MIMO H K turbulent channels as 
a function of the average SNR, fi, for various number of receiving apertures, 
N. Simulation Parameters: A] , n = 2, A 2 : n = 1, 0i,n = ?r/3, 02, n = tt/4, 

ckti — 2, 6o,n — 2. 



Fig. 2. ABEP of uncoded OSM for 2 X 2 and 2 X 4 MIMO H 
K turbulent channels as a function of the average SNR, //., for various 
values of link distances, L. Simulation Parameters: A = 1550nm, C 2 = 
1.7 X 10 -14 m -2 / 3 ,0i,„ = ir/3, 02,n = 7r/4. 


of the previously mentioned expressions, comparisons with 
complementary Monte Carlo simulated performance results 
are also included in these figures. As far as the performance 
of coded OSM systems is concerned, ABEP upper bounds vs. 
SNR have been obtained using (f43l) with (127b (see Fig. 0). 

Rg. m presents the ABEP performance as a function of 
the average SNR, /j,, of 2 x N MIMO OSM systems with 
N £ {1,2,3,4}. Independent and identically distributed 
branches are considered with Ai, n = 2, A 2 , ra = 1, $i, n = 7 r/ 3 , 
6*2 , n = 7 t/4 , a n = 2, 6 o, ra = 2. The obtained results clearly 
indicate that the ABEP curves, obtained using (l26l >. are in close 



Average SNR, p, (dB) 

Fig. 3. ABEP Comparison of 2 X 2 OSM with 1x2 coherent MRC systems 
employing DPSK, as a function of the average SNR, //, for various values of 
Ai n . Simulation Parameters: A 2 ,„ = 0, 0i, n = 0, 02 , n = 0. a n = 1.5, 
bo ,n — 1.5. 



Fig. 4. ABEP upper bounds of convolutional coded OSM for 2x2 and 
2x1 H-K turbulent channels as a function of the average SNR, //. for various 
values of link distances, L. Simulation Parameters: A = 1550nm, C 2 = 
1.7 X 10 -14 m -2 / 3 , 0i , n = 7r/3, 02 ,n = 7t/4. 


agreement with those obtained via simulations, verifying the 
correctness of the proposed analysis. Moreover, it is evident 
that the asymptotic ABEP curves correctly predict the diversity 
gain of the considered system for all tested cases. 

In Fig. [2] the dependence on the link distance L of the 
ABEP of a 2 x N MIMO OSM system is illustrated. The 
considered system is again equipped with either N = 2 or 
N = 4 receiving apertures and identically distributed branches 
are assumed. The parameters of the H-K distribution are 
calculated from (fl5] > and (fl~6]> assuming spherical wave propa¬ 
gation. Following |39l , it is further assumed that the operating 
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wavelength is A = 1550 nm and C 2 = 1.7 x 10~ 14 m _2 / 3 . As 
expected, the error performance deteriorates as L increases 
from L = 500m to L = 1500m. Moreover, it is evident 
that an increase in L from 500m to 1000m results in a more 
severe performance deterioration than in the case where L 
increases from 1000m to 1500m. In all cases considered, the 
analytical results obtained using ( l26l > are compared with the 
equivalent results obtained by means of Monte-Carlo computer 
simulations and again match very well. 

Next we compare the proposed OSM system with two alter¬ 
native coherent FSO systems that can provide performance en¬ 
hancements by means of transmit (MISO) or receive diversity 
(SIMO). It is noted that for similar aperture configurations, a 
fair comparison between coherent and IM/DD systems seems 
difficult as the same received laser power leads to different 
SNRs for each of these schemes ll23l . On the other hand, in 
order to perform a fair comparison between OSM and the 
alternative MISO or SIMO systems under the same propa¬ 
gation channel conditions, the aperture configuration of the 
FSO systems under comparison should be selected carefully. 
Specifically, because of the fact that the diversity gain of OSM 
equals to only the number of the receive apertures only, i.e. 
no transmit diversity gain is provided, the number of transmit 
or receive apertures of the alternative systems must be hence 
selected to be equal to the number of receive apertures of the 
OSM system. To this end, for a fair comparison in our paper 
a 2 x 2 OSM system is compared with the following two 
alternative FSO communication systems which also employ 
coherent detection: i) A 1 x 2 heterodyne FSO communication 
system which employs Differential Phase Shift Keying 1251 
and MRC or SC. 

ii) A 2 x 1 coherent FSO system employing the Alamouti 
scheme li22l and Binary Phase Shift Keying (BPSK). 

The instantaneous SNR at the output of the coherent MRC 
receiver assuming equal average SNR per receiving aperture, 
fi can be expressed as |29j 

N 

7MRC = In (46) 

n= 1 

whereas for SC is 

7sc = max{/i/i,/i/ 2 }. (47) 


In order to simplify the underlying mathematical analysis, it is 
assumed that the PDF of I n is given by © with the parameters 
A n being all zero, i.e. the PDF is the K-distribution. Thus, 
A4 r n (ji) can be readily obtained in closed form from (l33l) 
by replacing bo, n with b( l n /2. In Fig. [3] the ABEP of 2 x 2 
MIMO OSM links is compared with the ABEP of 1 x 2 
coherent FSO systems with DPSK. and identically distributed 
links are considered. In order to compare these systems under 
the same propagation conditions, it is assumed that a n = 1, 5, 
bo,n = 1.5, A 2 iU = 0 and A= {0,1,2,3}. As it can be 
observed, when either MRC or SC are employed, although 
coherent DPSK performs worse than OSM for values of Ai >n 
up to approximately 1, it outperforms OSM at lower values of 
A\ tn . Moreover, although the OSM outperforms the Alamouti 
scheme for Ai >n = 2 and 3, it performs similarly for high 
SNR values when A± tn = 1. It is noted that for Ai tU = 1 and 
lower values of A\ tn the Alamouti scheme yields the best 
performance of the considered OSM schemes. When more 
transmit appertures are employed, however, this advantage is 
compensated by the superior spectral efficiency of OSM and 
its lower hardware complexity as compared to coherent MRC. 
Specifically, as pointed out in m, OSM offers increased 
spectral efficiency by a factor log 2 (M). Moreover, as only one 
transmitting aperture is activated at any symbol duration, OSM 
has a lower decoding complexity as compared to conventional 
MRC and Alamouti schemes. 

In Fig. [0 upper bounds on the ABEP of convolutional 
coded 2x1 and 2x1 OSM systems are depicted, assuming 
similar propagation conditions to those considered in Fig. [2] 
Considering a convolutional code with rate 1/3 and constraint 
length of 3, its transfer function is given as Hi Eq. (8.2.6)] 


T[D(0),N] 


D(9fN 

l-2ND(9) 2 ' 


(50) 


Substituting ( l50l ) to ( |43| ). a union bound on the ABEP can be 
obtained as 


fir/2 


Pub < 


D(6») e 


TT\og 2 (M)J 0 (1-2D(9) 2 ) 2 


d 9. 


(51) 


The performance results of Fig. [0 clearly show that, as ex¬ 
pected, the incorporation of convolutional coding significantly 
enhances the performance of OSM systems, even when a small 
number of receive apertures is employed, even for N = 1. 


For MRC case, the ABEP can be deduced as |29l 

1 N 

p e = - 2 U^M- (48) 

n= 1 

For SC case, an analytical expression for the ABEP is more 
difficult to be deduced and, therefore, ABEP will be evaluated 
by means of Monte Carlo simulation only. 

As far as the Alamouti scheme is concerned, the instan¬ 
taneous SNR at the input of the demodulator of the optical 
receiver has a similar form as (l46l) |22| . For this scheme, the 
ABEP of BPSK can be evaluated as 



VI. Conclusion 

In this paper, the use of spatial modulation technique for 
coherent FSO communication systems has been proposed. 
We have provided a comprehensive analytical framework for 
error performance analysis in the presence of atmospheric 
turbulence scattering channel models which include the H-K 
distribution. The proposed framework reveals important in¬ 
formation about the performance of OSM over such turbulent 
channels, including the effect of fading severity and the achiev¬ 
able diversity gain. It also provides valuable insight into the 
impact of channel parameters on performance of OSM. Upper 
bounds for the ABEP performance of coded OSM systems 
have also been derived, demonstrating that coding techniques 
can greatly enhance the performance of OSM. Extensive 
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computer simulation performance evaluation results have been 
also obtained which have verified the accuracy of the analytical 
approach. Important trends about the performance of OSM 
for a variety of atmospheric turbulent scenarios and MIMO 
setups have also been identified. For example, it was shown 
that OSM can provide significant performance enhancements 
in the presence of atmospheric turbulence. The improvements 
are comparable to the ones offered by conventional coherent 
systems with spatial diversity, while outperforming the latter in 
terms of spectral efficiency and hardware complexity. Besides, 
under specific propagation conditions, OSM can yield bet¬ 
ter performance than conventional SIMO systems employing 
MRC or SC. We believe that the proposed framework is a 
useful tool for understanding the performance trend, important 
properties and tradeoffs of outdoor OSM operating in the 
presence of atmospheric turbulence. 
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